Algorithms for Enhanced Inter Cell Interference 
Coordination (elCIC) in LTE HetNets 

Supratim Deb, Pantelis Monogioudis, Jerzy Miemik, James P. Seymour 



Abstract — The success of LTE Heterogeneous Networks (Het- 
Nets) with macro cells and pico cells critically depends on 
efficient spectrum sharing between high-power macros and low- 
power picos. Two important challenges in this context are, (i) 
determining the amount of radio resources that macro cells 
should offer to pico cells, and (ii) determining the association 
rules that decide which UEs should associate with picos. In 
this paper, we develop a novel algorithm to solve these two 
coupled problems in a joint manner. Our algorithm has provable 
guarantee, and furthermore, it accounts for network topology, 
traffic load, and macro-pico interference map. Our solution is 
standard compliant and can be implemented using the notion of 
Almost Blank Subframes (ABS) and Cell Selection Bias (CSB) 
proposed by LTE standards. We also show extensive evaluations 
using RF plan from a real network and discuss SON based elCIC 
implementation. 

Index Terms — 4G LTE, Heterogeneous Cellular Systems, 
elCIC, Self-Optimized Networking (SON) 

L Introduction 

Wireless data traffic has seen proHfic growth in recent years 
due to new generation of wireless gadgets (e.g., smartphones, 
tablets, machine-to-machine communications) and also due to 
fundamental shift in traffic pattern from being data-centric 
to video-centric. Addressing this rapid growth in wireless 
data calls for making available radio spectrum as spectrally- 
efficient as possible. A key centerpiece is making the radio 
spectrum efficient is LTE heterogeneous networks (LTE Het- 
Net) or small cell networks |21|. In a HetNet architecture, in 
addition to usual macro cells, wireless access is also provided 
through low-powered and low-cost radio access nodes that 
have coverage radius around 10 m-300 m|7|. Small cells in 
LTE networks is a general term used to refer to Femto cells 
and Pico cells. Femto cells are typically for indoor use with a 
coverage radius of few tens of meters and its use is restricted to 
a handful of users in closed subscriber group. Pico cells have 
a coverage of couple of hundreds of meters and pico cells 
are open subscriber group cells with access permission to all 
subscribers of the operator. Picos are typically deployed near 
malls, offices, business localities with dense mobile usage etc. 
Picos are mostly deployed outdoor but there could be indoor 
deployments in large establishments etc. However, in LTE, 
since pico cells typically share the frequency band as macro 
cells, the performance of a low-power pico access node could 
be severely impaired by interference from a high power macro 
access node. Addressing this interference management riddle 
is key to realize the true potential of a LTE HetNet deployment 
and is the goal of this work. This work focuses on resource 
sharing between macro cells and pico cells. Note that macros 
and picos are both deployed in a planned manner by cellular 
operators. 
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Fig. 1. A typical LTE HetNet architecture with Macro and Pico access 
nodes. Pico-1 is used for throughput enhancement in a possible traffic hotspot 
location, Pico-2 and Pico-3 are used for improving edge throughput. 

A typical HetNet with pico and macro access nodes is 
shown in Figure [T] The high-power macro network nodes 
are deployed for blanket coverage of urban, suburban, or 
rural areas; whereas, the pico nodes with small RF coverage 
areas aim to complement the macro network nodes for filling 
coverage holes or enhancing throughput. There are two factors 
that could handicap the net capacity of a pico access node in 
the downlink. Firstly, the downlink pico transmissions to its 
associated UEs could be severely interfered by high power 
macro transmissions [llj. For e.g., in Figure [T] downlink 
transmissions to UEs associated with Pico-1 could easily be 
interfered by downlink transmissions of Macro- 1. Secondly, 
UEs, who are close to pico and could benefit from associating 
with a pico access node, could actually end up associating with 
the macro access node due to higher received signal strength 
from the high power macro access node[^ For e.g., UEs not 
too close to Pico-3 but still within the coverage area of Pico- 
3 could end up associating with Macro-2 because of higher 
received signal strength from Macro-2. Indeed, this could 
leave the pico underutilized and thus defeating the purpose of 
deploying that pico. Note that, it is the downlink interference at 
the pico UEs that needs additional protection from the macros; 
the uplink interference at the picos can be mitigated using the 
same power control principle in a macro only LTE network 161. 
Thus, for a pico cell based HetNet deployment to realize the 
promised theoretical gains, there are two important questions 
that need to be answered: 

1) How should downlink radio resources be shared so that 
pico UEs are guaranteed a fair share of throughput? 
Clearly, one needs to ensure that the pico transmissions 
are not badly hit by interference from macros. 

2) How to decide which UEs get associated with picos? 
Clearly, association based on highest signal strength is 
inadequate to address this challenge. 

This paper provides answers to these two coupled questions. 
Realizing the need to protect downlink pico transmissions by 
mitigating interference from neighboring macro cells, 3 GPP 
has proposed the notion of enhanced inter cell interference 

^In LTE networks, UEs associate typically with the cell with highest 
Received Signal Reference Power (RSRP). RSRP is a measure of the received 
signal strength of a cell at a UE and it is measured based on the strength of 
certain reference signals that cells broadcast. 
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coordination (elCIC) that provide means for macro and pico 
access nodes to time-share the radio resources for downlink 
transmissions. In simple terms elCIC standards propose two 
techniques. Firstly, each macro remains silent for certain 
periods, termed Almost Blank Subframes (ABS periods), over 
which pico can transmit at reduced interference. Secondly, 
the received signal strength based UE association in LTE is 
allowed to be biased towards the pico by a suitable margin. The 
details of how to set these ABS periods and how much to bias 
the association in favor of picos are left unspecified. In this 
paper, we answer these questions. We design our algorithms 
to meet the following goals: network- wide high performance, 
adaptability to network settings like propagation map and 
network topology etc., and scalability. 

A. Our Contributions 

In this work, we make the following contributions: 
1) Framework for network dependent elCIC: To the best of 
our knowledge, ours is the first work to provide a for- 
mal framework for optimizing Almost Blank Subframes 
(ABS) and UE-association in every cell by accounting 
for cell specific UE (load) locations, propagation map 
of each cell, macr-pico interference maps, and network 
topology. We also establish that computing the optimal 
solution with respect to maximizing a network utility is 
computationally hard. 

Efficient elCIC Algorithms: We next provide an efficient 
algorithm to compute ABS and UE- associations (and 
corresponding CSB) in an LTE HetNet. Our algorithm 
is provably within a constant factor of the optimal and 
scales linearly with the number of cells. Furthermore, 
our algorithm is amenable to distributed implementation. 
Evaluation using Real RE Plan: We perform extensive 
evaluation of our algorithm on a Radio-Frequency map 
from a real LTE deployment in New York City and 
demonstrate the gains. The results show that, our al- 
gorithm performs within 90% of the optimal for real- 
istic deployment scenarios, and, 5^^ percentile of UE 
throughput in the pico coverage area can improve up 
to more than 50% compared to no elCIC; the improve- 
ments can be 2x for lower throughput percentiles. 
Practical Eeasibility with SON: Finally, we discuss 
the challenges of implementing elCIC within Self- 
Optimizing (SON) framework and describe a prototype 
along with the associated challenges. 
The rest of the paper is organized as follows. Section |ll| 
provides a background on elCIC and describes some impor- 
tant related work. In Section |lll| we describe our network 
model. Section |X] states the problem and formally derives 
the computational limits of the problem. Our main algorithm 
for jointly optimizing ABS parameters for each cell and 
macro/pico association for each UE is provided in Section [V]- 
VII Section VIII describes how a given choice of UE associ- 
ation can be translated into cell selection bias parameters and 
also how ABS numbers can be converted into ABS patterns. 
Section IX presents evaluations using RF plan from a real LTE 
deployment. Finally, Section |X| discusses SON based elCIC 
implementation along with a prototype. 

II. Background: eICIC and Related Work 
A. Enhanced Inter Cell Interference Coordination (elCIC) 

The elCIC proposal in LTE standards serves two important 
purposes: allow for time-sharing of spectrum resources (for 
downlink transmissions) between macros and picos so as 
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Fig. 2. An illustration of how an LTE frame can consist of ABS subframes. 
A pico can transmit over ABS subframes with very little interference from 
macro, and it can also transmit over any other non-ABS subframe when it 
receives high interference from macro. 



to mitigate interference to pico in the downlink, and, allow 
for flexibility in UE association so that picos are neither 
underutilized nor overloaded. In elCIC, a macro eNodeB can 
inject silence periods in its transmission schedule from time 
to time, so that interfering pico eNodeB s can use those silence 
periods for downlink transmissions. Furthermore, to ensure 
that sufficient number of UEs get associated with a pico, 
the elCIC mechanism allow UEs to bias its association to 
a pico. Before we discuss these mechanisms in more details, 
we provide a very brief introduction to format of downlink 
transmissions in LTE. 

Downlink transmission format: In LTE, transmissions are 
scheduled once every subframe of duration 1 ms; 10 such 
subframes consist of a frame of length 10 ms. Each subframe 
is further divided into 2 slots of duration 0.5 ms each. Each 
slot consists of 7 OFDM A symbols. While we do not need 
any further details for our discussion, the interested reader 
can refer to |17| for extensive details of LTE downlink 
transmission format. We now describe two important features 
of elCIC. 

Almost Blank Subframes (ABS): In order to assist pico 
downlink transmissions, the macro eNodeB s can mute all 
downlink transmissions to its UEs in certain subframes termed 
almost blank subframes (ABS). These subframes are called 
"almost blank" because a macro can still transmit some 
broadcast signals over these subframes. Since these broadcast 
signals only occupy a small fraction of the OFDMA sub- 
carriers, the overall interference a macro causes to a pico 
is much less during these ABS periods. Thus, the pico can 
transmit to its UEs at a much higher data rate during ABS 
periods. Note that, a pico is also allowed to transmit to its 
UEs during non-ABS periods. This could provide good enough 
performance to UEs very close to the pico. An example of 
ABS schedule is shown in Figure [2] 

Flexible User Association and Cell selection bias (CSB): 
Typically in cellular networks, when a UE device (UE) has 
to select a suitable cell for association, it chooses the one 
with maximum received signal strength. However, if the same 
strategy is extended to HetNet deployments with both macro 
and pico cells, this could lead to underutilization of the 
pico eNodeB 's. This is because, picos transmit at very low 
power and thus, unless a UE is very close to the pico, signal 
strength from the macro is likely to be larger for the UE. 
To overcome this, LTE standards have proposed a concept 
called cell selection bias which works in the following manner. 
Suppose the cell selection bias of cell-i is a^. Denote by Pi 
as the reference-signal received power (in dBm) from cell-z 
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as measured by a UE. Here a cell could be a pico-cell or a 
macro-cell. Then, the UE associates with cell-/c such that 

k = argmax(Pi + a^). 

i 

Thus, by assigning larger (smaller) bias to picos compared to 
macros, once can ensure that the picos are not underutilized 
(over-utilized). The bias values are broadcast by the cells to 
assist UEs make the right association decision. 

B. Connection with Other Interference Mitigation Techniques 

It is instructive to discuss the connection of elCIC with two 
other interference mitigation techniques found in LTE specifi- 
cations, namely, the frequency domain Inter-cell Interference 
Coordination (ICIC) and Coordinated Multi-Point (CoMP). 

Connection with ICIC: ICIC pre-existed elCIC in terms 
of both standards development and initial deployments. In 
ICIC based schemes, cell resources are divided into fre- 
quency bands and transmit power profiles to reduce inter- 
cell interference |16|, |18|, |9|. The obvious difference is 
that ICIC works in the power spectral density (PSD) domain 
as opposed to the time domain of elCIC. ICIC has certain 
limitations compared to elCIC. In ICIC based schemes, the 
different power profiles can either span different sub-bands 
(thus leading to a non-uniform frequency re-use factor) or 
it can span multiple carriers. The disadvantage of sub-band 
ICIC is that standards did not include the capability of varying 
the reference signal power according to the sub-band, giving 
rise to data demodulation issues for QAM-modulated symbols. 
Another challenge in deploying sub-band based ICIC is to 
find sufficient number of PSD patterns in a dense deployment 
with macros and picos. The carrier based ICIC, on the other 
hand, introduces significant number of inter-carrier handovers. 
The requirement to define multiple carriers also increases the 
cost of small cells. Due to these reasons, elCIC has attracted 
much attention of operators and standard bodies especially for 
mitigating interference to picos in HetNets. Having said that, 
ICIC based macros and elCIC are not either-or propositions 
and elCIC benefits can come on the top of existing ICIC 
based macro deployments. Our model and framework easily 
incorporates existing ICIC based macros; see Remark |2] in 
Section |7^/6>r details. 

Connection with CoMP: Coordinated Multipoint Access 
(CoMP) is also another edge-rate improving technology for 
downlink and uplink both. In downlink, CoMP can be 
considered an extension of multiuser MIMO (MU-MIMO), 
that achieves interference mitigation by transmitting simul- 
taneously from multiple cells with properly chosen antenna 
weights so as to achieve some optimal physical-layer oriented 
metric like zero-forcing or joint Minimum Mean Square Error 
(MMSE). In spite of theoretical performance benefits, CoMP 
faces severe practical challenges due to high backhaul re- 
quirements between base stations, extremely tight time-phase- 
frequency synchronization across the set of collaborating cells, 
and requirement of pico cells to be equipped with antenna 
arrays leading to cost increase of supposedly low-cost picos. 
Evidently, it will take significant investment and time for 
CoMP to be widely deployed in LTE networks that are still in 
their early days. elCIC on the other hand is much simpler to 
deploy and is likely to be adopted much sooner. To this end, 
our work optimizes elCIC which does not account for CoMP. 
Nevertheless, jointly optimizing a modified version of elCIC 
to account for CoMP is a technically challenging problem and 
is left as a future work. See Remark |3] in Section [Iv] for further 
details. 



C. Related Work 

The elCIC proposal is relatively new for LTE Hetero- 
geneous networks. In |7|, the authors present a very good 
introduction to the concept of elCIC in LTE HetNets. In ifTTIl . 
the authors provide an excellent survey on elCIC and the 
motivation behind elCIC proposal in LTE standards. In a 
recent work |12|, the authors present simulation studies to 
understand the dependence between network performance and 
elCIC parameters. However, the authors primarily consider 
uniform elCIC parameter in all the cells; clearly, the right 
choice of elCIC parameters should vary across cells and 
account for propagation map, cell-load etc. Also, the authors 
do not present a framework to optimize the elCIC parameters. 

In the previous subsection, we have talked about another 
interference mitigation technique called ICIC that has been an 
area of active research in the recent past. In 1 16], 1 18 1, |9 |, the 
impact of ICIC has been studied for LTE and LTE HetNets. 
The concept of soft-frequency reuse has been formally studied 
in (201, where, the authors optimize downlink transmit power 
profiles in different frequency bands. The work closest to 
ours in principle is |13 | that considers the problem of UE 
association and ICIC in a joint manner. 

III. System Model 

A. Terminologies 

In addition to the notion of subframes described in Sec- 
tion [ll| we will use the following terminologies in this paper. 

UE (user equipment): UE refers to the mobile device. 

eNodeB (eNB): The eNB|^ is the network element that 
interfaces with the UE and it performs radio resource manage- 
ment, admission control, scheduling, QoS enforcement, cell 
information broadcast etc. It hosts critical protocol layers like 
PHY, MAC, and Radio Link Control (RLC) etc. 

Macro cell: In LTE heterogeneous networks (HetNet), a 
macro cell has a base station transmitter with high transmission 
power (typically 20 W-40 W), sectorized directional antennas, 
high elevation, and thus ensuring a cell coverage radius 
typically around 0.5 km-2 km. 

Pico cell: As opposed to a macro cell, a pico transmitter 
is characterized my much lower transmission power (typically 
2 W-5 W), omnidirectional antennas, low antenna height, low 
cost, and has a cell coverage radius of around 100-300 m. 
Pico cells are underlayed on the macro-cellular network to 
fill coverage holes and to enhance capacity in traffic hotspot 
locations. 

Reference Signal Received Power (RSRP): Every UE in LTE 
makes certain measurements of received signal strength of all 
nearby cell transmitters. RSRP is the average received power 
of all downlink reference signals across the entire bandwidth 
as measured by a UE. RSRP is taken as a measure of the 
received signal strength of a cell transmitter at a UE. 

B. Network Model and Interference graph 

Since the elCIC proposal by LTE standard aims to protect 
downlink pico transmissionqj and our goal is to develop 
solutions for optimal elCIC setting, we only consider downlink 
transmissions in this work. 

Network Topology: Our system model consists of a network 
of macro and pico (also called pico in this paper) eNBs. M 

^eNB is equivalent to base station in traditional cellular voice networks but 
it has more functionalities. 

^The uplink problem (via power control |6|) in presence of picos is not 
different from macro only network because UE capabilities from a transmit 
power point of view remain same. 
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TABLE I 

List of Parameters and Key Optimization Variables 



be 



Notation 


Description 


U, u, N 

{u e u) 


Set of UEs, index for a typical 
UE, number of UEs, respectively 


M, m, M 
(m G M) 


Set of macros, index for a typical 
macro, number of macros, respectively 


1 1 Pi ^ 


Set of picos, index for a typical 
pico, number of picos, respectively 


ruu 


The macro that is best for UE-w 


' u 


Data-rate achievable by UE-ii from JTiu 
(in bits/sub-frame) 




The pico that is best for UE-ii 


pico^ABS 


Data-rate achievable by UE-w from pu when all 
interrering macros are muted 
(in bits/sub-frame) 




Data-rate achievable by UE-w from pu when 
all/some interfering macros are transmitting 
(in bits/sub-frame) 




Set of macro eNB's that interfere 
with pico p 


lAm 


Set or UEs tor whom macro-m is the 
best macro eNB 


7 J 

Up 


Set of UEs for whom pico-^ is the 
best pico eNB 


Ap 


Variable for ABS subframes used/received by pico-p 


Nm 


Variable non-ABS subframes used by macrio-m 




Variable denoting UE-ii's air- time from macro 




Variables denoting UE-w's air- time from pico 
over ABS and non-ABS subframes respectively 


Ru 


Variable denote UE-w's average throughput 


z, p 


Vector of all primal and dual variables, respectively 



denotes the set of macros and V denotes the set of picos. We 
also use m and p to denote a typical macro and a typical pico 
respectively. 

Interference Modeling and Macro-pico Interference Graph: 
We now describe our interference model. For the purpose of 
elCIC algorithms, it is important to distinguish macro-pico 
interference from the rest. 

• Macro-pico interference: For each pico-p, the set of 
macros that interfere with it is denoted by 2p C . The 
macros in the set 2p need to be silent during any ABS 
subframes used by pico-p. Thus, UEs of pico-p can be 
interfered by m G 2p only during non-ABS subframes. 

• Macro-macro and pico-pico interference: Due to 1 : 1 fre- 
quency re-use in most LTE networks, picos can interfere 
with each other and similarly for the macros. A pico UE 
can be interfered by another pico both during ABS and 
non-ABS subframes. 

To better understand the distinction between the two kinds 
of interferences from an elCIC point of view, consider a 
pico associated UE-i^'s interfering cells. Suppose the total 
interference power it receives from all other interfering picos 
and all interfering macros (those in the set Xp) be PfnTi^) 
and PY^^^^{u), respectively. Denoting by Prx{u) the received 
downlink power of UE-i^, the downlink SINR of UE u E Up 
can be modeled as 



S\m{u) = 



(1) 



for ABS sub-frames 
for non-ABS sub-frames 



This is because, during ABS sub-frames all interfering macros 
of pico-p remain silent and so the only interference is from the 
interfering picos of p. However, during non-ABS sub-frames, 
there is interference from all interfering picos and macros both. 
Instead, if UE-ii were a macro-UE, the SINR expression would 



SINR(ii) 



Prx{u) 



(2) 



where Pf^Ti^) Pj^T^^i^) denote the interference from 
interfering picos and macros respectively. 

Thus the interference graph relevant from picos ' point of 
view is the bipartite graph formed by joining edges from any 
p to macros in the set Ip. The graph neighbors of a pico-p 
should all remain silent during ABS periods usable by picos. 

Remark 1. The elements oflp can be obtained either through 
cell- adjacency relationship or based on whether the received 
signal from macro eNB to p is above a threshold. 

User model: To start with, we will consider a scenario 
where there is a set of static UE's denoted by U, and also, 
we know for each UE-u ( i) the best candidate macro in terms 
of RSRP and average PHY data-rate from the macro r^^^^^ , 
(ii) the best candidate pico, if any, and average PHY data- 
rate in ABS and non-ABS subframes given by 7-^*^0, abs" 
and rJJ*^^ respectively. Note that, the value of r^^^^^^ can 
be obtained from the SINR expression ^ using LTE table 
lookup for conversion from SINR to rate, and similarly the 
values of rJJ^^^'^^*^ and r^^'^'' from the SINR expression ([T]). 
Alternatively, one can use Shannon capacity formula (with 
some offset in SINR) to obtain r'^acro ^^pico,ABS ^^pico f^^^ 

the corresponding SINR expressions. Clearly, the average PHY 
data rate that a UE receives from pico is higher in ABS frames 
due to reduced interference from nearby macros. In fact, the 
average PHY data rate from a pico in a non-ABS subframe 
is likely to be very small in many instances due to very high 
interference from macro. 

Note that picos could be deployed indoor or outdoor (typi- 
cally outdoor) but makes no difference to our framework. Only 
the pico to UE propagation (and hence data rates) change ap- 
propriately. The main parameters and some of the optimization 
variables (introduced later) are captured in Table |l] 

IV. Problem Statement and Computational 
Hardness 

We will first develop an algorithm to find optimal ABS and 
CSB configuration with static UEs scenario where we have 
the precise knowledge about number of UEs in different cells 
along with PHY data rates. We will describe in Section |X| 
how Monte -Carlo based techniques can be used along with 
our algorithm for this scenario where UE densities and SINR 
distributions are known instead of exact UE locations. 

The essence of elCIC approach is to compute optimal 
association (either to macro or to a pico) rules for UEs, and 
also compute how macros and picos share radio resources in 
time domain. Thus, we will first formulate a problem for the 
optimal choice of (i) UE association, i.e., which UEs associate 
with the best macro and which ones associate with the best 
pico, (ii) the number of ABS subframes reserved for interfered 
picos by each macro eNB. We will denote by Nsf as the total 
number of subframes over which ABS subframes are reserved 
(typically, Ngf = 40). We will also refer to the quantity Ngf 
as ABS-period. 

Optimization variables: Let Nm be the number of sub- 
frames for which macro-m can transmit during each ABS- 
period (clearly, Nsf — N^, ABS subframes are offered by 
macro-m in each ABS-period). Let Xu be the time-average 
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air-timqjin sub-frames per ABS-period UE-u gets from rriu, 
the candidate best macro for UE-u. Note that Xu need not be an 
integer. The airtime that UE-u gets from a pico can be during 
ABS subframes or regular subframes because pico eNBs 
can transmit during ABS subframes and regular subframes. 
To this end, we define and y'^^ as the time-average 
airtime in subframes per ABS-period UE-u gets from pico Pu 
(best candidate pico of UE-u) in ABS subframes and regular 
subframes respectively. Also, let Ru be the average throughput 
UE-u achieves. 

Optimization constraints: These are explained below. 

1) Association constraints: The association constraint es- 
sentially states that a UE can associate with either the 
macro or a pico but not both. Thus, 

\/ueU : XuiVu + Vu"^) = (3) 
Xu >0,yu> 0, yu"^ > . 

so that, either the total airtime u gets from macro is 
zero or the total airtime u gets from pico is zero. 

2) Throughput constraints: This implies that the average 
throughput for a UE-u, Ru, cannot be more than what 
is available based on the air-times from associated 
macro/pico. 

w ^ 7 / r> ^ macro . pico.ABS A . pico nA 

v u eu : Ru < Xu + K yu +K yu • 

(4) 

Note that, from ([3j, if UE-u is associated to a macro, 
then Ru < r^^^^^x^; if UE-u is associated to a pico, 
then Ru < rJJ^^^'^^*^^;^ + r^^^^T/^J^. 

3) Interference constraints: The interference constraint 
states that the ABS subframes used by a pico p are 
offered by all macros in the set 2p that interfere with 
the pico. In other words, 

V (p, m G Xp) : Ap + N^< Nsf . (5) 



4) Total airtime constraints: This ensures that the total 
time-average airtime allocated to UE's from a macro or 
a pico is less than the total usable subframes. This can 
be described using the following inequalities. 

V m G M : ^ X. < AT^ (6) 

ueUm 

V p G P : ^ 2/^ < (7) 

u^Up 

VpGP : ^ (y^ + 2/:-*) < iV,/ (8) 

u^Up 

Optimization objective and problem statement: The opti- 
mization objective we choose is the tried and tested weighted 
proportional-fair objective which maximizes ^^Wul^Ru, 
where Wu represents a weight associated with UE u. This 
choice of objective has three benefits. Firstly, it is well known 
that proportional-fair objective strikes a very good balance 
between system throughput and UE-throughput fairness (221 . 
This bodes well with the goal of improving cell-edge through- 
put using elCIC. Secondly, such a choice of objective gels very 
well with the underlying LTE MAC where the most prevalent 
approach is to maximize a proportional-fair metric. Finally, 
the weights Wu provides a means for service-differentiation [2 1 
which is a key element in LTE. The weights may be induced 

^This time-average airtime can be achieved through MAC scheduling, 
particularly weighted proportional-fair scheduling. We are implicitly assuming 
two time-scales here: the ABS selection happens at a slow time-scale and 
MAC scheduling happens at a fast time-scale resulting in a time-average 
airtime for each UE. 



from policy-determining functions within or outside of radio 
access network. Though we use proportional-fair metric in this 
paper, our algorithms can be easily modified to work for other 
utility functions. 

The problem can be stated as follows: 



We also call the optimization objective as system utility and 
denote it, as a function of the all UE's throughput- vector R, 
by 

Util(R) = ^ w;„ In i?„ . 

U 

Remark 2. f Accounting for ICICj Though elCIC or time- 
domain resource sharing is the preferred mode of resource 
sharing between ma cro a nd pico cells in LTE for reasons 
mentioned in Section \ITB\ elCIC could co- exist with ICIC in 
macro-cells. Our framework can easily account for this with 
only changes in the input to the problem. In ICIC, the OFDMA 
sub-carriers of a macro are partitioned into two parts: low- 
power sub-carriers of total bandwidth Bi and high-power sub- 
carriers of total bandwidth B^. Thus for a UE u that receives 
signal from macro-m, the spectral efficiency over low-power 
subcarriers (say, rfu'^) is different from spectral efficiency over 
high-power subcarriers (say, r]^^^^). This causes the downlink 
rate from the candidate parent macro, r JJ^^^^^, to be expressed 
as r;;^^^^^ = r]^'^^Bh + rfu"^ Bi. Furthermore, the reduced 
interference from ICIC-using macros changes the macro-pico 
interference graph structure and pico to UE non-ABS rates. 
Importantly, only the input to OPT-ABS has to be modified 
to account for ICIC in frequency domain. 

Remark 3. fCoMPj Our framework assumes no CoMP based 
deployments. This would be the case in the most LTE deploy- 
ments in the foreseeable futu re du e to practical challenges 
of CoMP outlined in Section \II-B\ Technically speaking, to 
optimize a modified version of elCIC that accounts for CoMP, 
the association constraint given by Q would not be required, 
the throughput constraint must account for collaborating 



OPT-ABS 

Given: A set of UEs a set of picos V, and a set of 
macros M. For each UE u eU are given the following: 
best candidate parent macro along with PHY data 
rate to rriu denoted by r^^^^^, the best candidate parent 
pico Pu along with ABS and non-ABS PHY data rate 
to candidate parent pico denoted by 7-^*^0,^55 ^pzco 
respectively. We are also given the macro-pico interference 
graph in which the interfering macros of pico-p is denoted 
by Xp. Finally, Ngj is the total number of subframes. 

To compute: We wish to compute the number of ABS 
subframes Ap each pico-p can use, the number of non-ABS 
subframes Nm left for macro-m 's usage, a binary decision 
on whether each UE-u associates with its candidate parent 
pico or candidate parent macro, throughput Ru of each UE- 
u, so that the following optimization problem is solved: 

maximizei^^ ,yA^yr.A^Ap,N^,R^}^uJu^n Ru 

u 

subject to, (|3|, 0, (|5|, (|6|, Q, (H 

V ip,m€lp) : Ap, iV™ eZ+ , 

where Z+ denotes the space of non-negative integers. 
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cells in CoMP, and the total airtime constraints must be 
modified to reflect CoMR We leave this as a future work. 

A. Computational hardness 

It can be shown that the ABS -optimization problem is NP- 
hard even with a single macro but multiple picos. We state the 
result as follows. 

Proposition 1. Even with a single pico and a single interfering 
macro, the OPT-ABS problem is NP-hard unless P = NP. 

Proof Follows by reducing the SUBSET-SUM problem (W\ to 
an instance of OPT-ABS problem. See|5|. □ 

In light of the above result, we can only hope to have algo- 
rithm that is provably a good approximation to the optimal. In 
the following, we will develop an algorithm with a constant- 
factor worst case guarantee; we show extensive simulation 
results to demonstrate that our algorithm is within 90% of 
the optimal is many practical scenarios of interest. 

V. Algorithm Overview 

Our approach to the problem is to solve it in two steps. 

1) Solving the relaxed NLP: In the first step, we solve 
the non-linear program (NLP) obtained by ignoring 
integrality constraints on Ap and N^n and also the 
constraint that a UE can receive data either from pico 
or macro but not both. Specifically, in this step, we 
maximize Util(R) subject to the constraints (|4])-([8]); and 
we allow Ap and Nm to take non-integer values. Note 
that ignoring the constraint ^ means that UEs can 
receive radio resources from macro and pico both. For 
notational convenience, we also denote the vector of 
constraints (|4])-([8]), in a compact form as gi?(.) < . 
Thus, the RELAXED- ABS problem can be denoted as, 

RELAXED - ABS : 

maximizei^^^^A^^^A^^^^^^I In Ru 

u 

subject to, gi?(.) ^ 

V [p.meXp) : Ap, G R+ , 

where R+ denotes the space of non-negative real num- 
bers and gi?(.) denotes the vector of constraints ([^-([s]). 

2) Integer rounding: In the second step, we appropriately 
round the output of the non-linear optimization to yield 
a solution to the original problem that is feasible. 

VI. Algorithm for Relaxed Non Linear Program 

Towards solving ABS -RELAXED, we use a dual based ap- 
proach |[T4l . Pi which has been successfully applied to many 
networking problems for its simplicity of implementation |4|, 
[19|. In the following, we show that, a dual based approach to 
our problem lends to a decomposition that greatly reduces the 
algorithmic complexity and also makes the approach amenable 
to distributed implementation while retaining the core essence. 

In a dual based approach, it is crucial to define a suitable 
notion of feasible sub- space so that the solution in each 
iteration is forced to lie within that sub-space. The choice of 
the sub- space also has implication on the convergence speed 
of the algorithm. To this end, we define the sub- space 11 as 
follows: 



n={x,y,A,N: Ar,j,Ar^<Ar,j, ^ Xu<N,f, 

ueUm 

We have used bold-face notations to denote vectors of 
variables. Clearly, any solution that satisfies the constraints 
described in the previous section lies within 11. In the follow- 
ing discussion, even without explicit mention, it is understood 
that optimization variables always lie in 11. 

We now describe the non-linear program (NLP) obtained 
by treating Ap and Nm as real numbers. The Lagrangian of 
the relaxed NLP can be expressed as follows: 

£(x, y, A, N, A, /it, /3, a) = ^ Wu In Ru (10) 

u 

E\ ( T? macro pico, ABS A pico nA\ 

u 

- ^ llp^m{Ap + Nm- Nsf) 
p,m^Xp 

P ueUp 

Additional notations: We use bold-face notations to express 
vectors. For example A denotes the vector of values A^^. 
The variables A, /i, /3, a's are dual variables and so called 
Lagrange-multipliers which also have a price interpretation. In 
the rest of the paper, p denotes the vector of all dual variables, 
i.e., p = (A, /i, /3, an Similarly, the variables x, y, A, N are 
referred to as primal variables and we use z to denote the 
vector of all primal variables, i.e., z = (x, y, A, N). Thus, we 
denote the Lagrangian by >C(z, p) and express it as 

/:(z,p)=Util(R) - p'gR{z) . 

The dual problem of RELAXED- ABS can be expressed as 

min7:>(p) , (11) 

p>0 

where, 

V{p) = max/:(z, p) . (12) 

zGll 

Since the RELAXED- ABS is a maximization problem with 
concave objective and convex feasible region, it follows that 
there is no duality gap |3|, and thus 

RELAXED-ABS Optimal = minP(p) . 

p 

Iterative steps: First the primal variables are initialized to 
any value with 11 and the dual variables are initialized to zero, 
and then, the following steps are iterated (we show the update 
for iteration-Q: 

1) Greedy primal update: The primal variables z^ in 
iteration- (t + 1) are set as 

zt+i = argmax£(z, p^) . (13) 

^The dual variables are often denoted by p because they have the interpre- 
tation of prices. 
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2) Subgradient descent based dual update: The dual vari- 
ables are updated in a gradient descent like manner as 



(14) 



where pt is the dual variable at iteration-t, 7 is the step- 
size, and [.]+ denotes component- wise projection into 
the space of non-negative real numbers. 
The above steps are continued for sufficiently large number 

of iterations T and the optimal solution to RELAXED-ABS 

is produced as 

T 

t=l 

The computation of greedy primal update is not immediate 
at a first glance, however, the subgradient descent based dual 
update step is straightforward. Thus, for the above algorithm 
to work, there are two important questions that need to be 
answered: (i) how can the greedy primal update step be 
performed efficiently? (ii) how should the step size 7 and the 
number of iterations T be chosen? In the following, we answer 
these questions. 

A. Greedy Primal Update: Decomposition Based Approach 
We now argue that the problem of computing 

argmax/:(z, pt_i) 

zGll 

can be decomposed into UE problem, macro problem and pico 
problem each of which is fairly straightforward. Towards this 
end, we rewrite £(z, p) as follows. 

u m 

p 

where, 

Fu{v^ Ru) = Wu \nRu - XuRu 



I / ^ yuK^uiu Hp LXp)^ / ^ Uu y^ui u ^p) 

u^lAp u^lAp 

It follows that, 
max£(z, p) 

= y^maxF^,(p,i^^,) + y^maxGm(p, {Xu}u^Urr^^ Nm) 
u m 

+ Xmaxi7p(p, {yu}ueUp,Ap) - Nsf 
p 

where the max in the above is with respect to appropriate 
primal variables from Xu^Vu^^ ^m's and Ap's. The above 
simplification shows that greedy primal update step can be 
broken up into sub-problems corresponding to individual UEs, 
individual macros, and individual picos; each of the sub- 
problems has a solution that is easy to compute as follows. 
We thus have the following. 



Greedy primal update: In iteration-t, the greedy primal 
updates are as follows: 

• User primal update: In greedy primal update step of 

iteration- for each UE-u, we maximize Fn(Pt, Ru) 

by choosing Ru{t -\- 1) as 



(15) 



Macro primal update: In greedy primal update step 
of iteration- (t + 1), for each macro-m, we maximize 

Gm{Pt, {^ujueUrr,, Nm) by choosing Nm as 

Nm{t + 1) = ^^/I{(/3^(^)-Ep:^eT^ /ip,^(t)>o)} • (16) 

To compute all {xu}ueUm^ ^^^h macro-m computes the 
best UE li* in iteration-t as 



arg max(A^(t)r^ 

ueUm 



Pm{t) > 0) 



where ties are broken at random. Macro-m then chooses 

Xu{t + G Um as 



Nsf for u = 
foru^ ix* 



(17) 



Pico primal update: In iteration- (t + 1), for each pico-p, 
we maximize Hp{pt, {VujueUpjAp) by choosing 



Nsfl 



5/M(/3p(t)-E^:^ex,Mp,^(t)>0)} • 



(18) 



To compute all {yu}ueu 's, each pico-p computes the 
current best UE u^iABS] and u^{nABS) as follows: 

ul{ABS) = argmax(A,rr^'^^^-/3p(t)-ap(t) > 0) , 



u^{nABS) — arg max(A^i(t) 



u^lAp 



ap(t)>{)) . 



where ties are broken at random. Pico-p then chooses 

yu{t + G Up as 



Nsf for u = ul{ABS) 

for u^ul{ABS) ^^^^ 



Similarly, we set y^^it) based on u^{nABS) as follows. 

nA(, , n _ / ^sf for u = u;{nABS) 

+ - \ for / u^inABS) ^^^^ 

B. Overall Algorithm for RELAXED-ABS 

We now summarize the algorithmic steps for solving 
RELAXED-ABS. Algorithm [T] formally describes our algo- 
rithm. 

We next derive the step-size and sufficient number of 
iterations in terms of the problem parameters. 

C. Step-size and Iteration Rule using Convergence Analysis 

Towards the goal of estimating the step- size and number of 
iterations, we adapt convergence analysis for a generic dual 
based algorithm is provided in 1 14 |. We show that the structure 
of ABS -RELAXED lends to a simple characterization of 
the step-size and number of iterations in terms of problem 
parameters. 

In this section, we denote by r^ax and rmin as the maxi- 
mum and minimum data rate of any UE respectively. We also 
denote by Wm as the total weight of all candidate UEs of 
macro-m and similarly for Wp. Also W denotes the vector of 
Wrn's and Wp's. Also Um, Up, Umax denote, respectively, the 
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Algorithm 1 Optimal RELAXED-ABS: Algorithm for 
Solving RELAXED-ABS 

1: Initialization: Initialize all the variables cc, i/, A, AT, A, /it, /3, a 

to any feasible value. 
2: for t = 0, 2, 3, . . . , T iterations do 

3: Primal update: Update the primal variables z (t) by using UE's 
update given by (T5] l, macro updates given by |T6| , (Tt] ), and 
pico updates givenby (Tsll, |T9l ( [20| . 

4: Z)wa/ Update: For each Ub-'i~Ati(T) is updated as 

Xu(t) ^ [Xu(t-1)+ 
j{Ru (t) -r~xu (t) -rr^^'^^^2/;f (t) -^^^^^^ (t) )] + . 
For each macro-m, we update its dual price /3m as follows: 

/3m(t) ^ [/3^(t - 1) + 7( ^ Xu{t) - Nm{t))]^ 

For each pico p, we update all dual variables /3p, ap and /ip,m 
for all m G Xp, as follows: 

(t) ^ K,m(t - 1) + 7(^p(i) + iVm(i) - Nsf)r 

Mt) ^ [Mt - 1) + 7( E vuit) - A,m+ 

a,{t) ^ [a,{t - 1) + 7( ^ {y^{t) + yl^{t)) - N,f)] + 
5: end for 

6: The optimal values of the NLP are obtained by averaging over 
all iterations: 

T 



number of candidate UEs in macro-m, number of candidate 
UEs in pico-p, and maximum number of UEs in any macro 
or pico. 

Proposition 2. Let T^t^T^t-,'^'' (VtiVt-, P*j denote the vector of 
primal (dual) variables at time t, averaged over all iterations 
fromO—t, and at optimality, respectively. Under mild technical 
assumptions JJ]/, we have the following: 

(i) D(pr)-I?(p')<^ + 2|! 



be used to set 7 and T as follows: 



where 



(ii) ra/(R*) - Util{tiT) < 



Q' = iV,^^(iVrL. + M + P + 2/) 



2 — (1 H~ '^^max + J^^^ ) 



Proof See [5|. 



□ 



Remark 4. (On the proof of Proposition [2j The main 
contribution of the proof of Proposition |2] is to show that the 
norm of optimal dual variable ||p* || can be upper bounded by 
network parameters. This upper bound, along with adaptation 
of convergence analysis in [14], readily characterizes the step- 
size and number of iterations for RELAXED-ABS; this is unlike 
arbitrary convex programs where the convergence results are 
in terms of a generic slater vector 47?!/ . 

Remark 5. (Step-size and Number of Iterations.) The 
step-size and number of iterations can be set based on the two 
following principles: 

1) Suppose we want the per-UE objective to deviate from 
the optimal by no more than e. Then, Proposition |2] can 



< 



Ne 



and 



27T 



< 



Ne 



The above imply. 
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Ne 



and T 



(QB. 
V Ne ) 



(21) 



(22) 



Since the maximum interferers Imax is typically a 
small number, a moments reflection shows that 7 = 
0{e/N^frl^ax) T = 0{wl^^^Umax/e^rmin) where 
u)max is the maximum vaue of Wu- In other words, the 
number of iterations simply depends on the maximum 
UEs in any cell and not on the overall number of UEs. 
2) The number of iterations required can be significantly 
reduced using the following observation. Suppose the 
macro-pico interference graph can be decomposed into 
several disjoint components. In that case, we can run 
the RELAXED-ABS algorithm for each component in- 
dependently (possibly parallaly). For each interference 
graph component, we can use the step-size and iteration 
rule prescribed in the previous paragraph. 

Nil. Integer Rounding of RELAXED-ABS 

In this section, we show how solution to RELAXED-ABS 
can be converted to a feasible solution for the original problem 
OPT-ABS. There are two challenges in performing this step. 
Firstly, in OPT-ABS, each UE can receive resources either 
from a macro or a pico but not both unlike RELAXED-ABS. 
Secondly, as with all dual based sub-gradient algorithms, after 
T iterations of running RELAXED-ABS, the solution may 
violate feasibility, albeit by a small margin |[T4l . Thus, we 
need to associate each UE with a macro or a pico and round 
the values of A/'^'s and Ap's so that the overall solution is 
feasible and has provable performance guarantee. 

To this end, we first introduce the following rounding 
function: 



RndAT^^ {x) 




X > 
X < 



(23) 



The rounding algorithm is formally described in Algorithm |2] 
The algorithm has three high-levels steps. In the UE as- 
sociation step, each UE who gets higher throughput from 
a macro in the solution of ABS-RELAXED is associated 
with a macro, and, each UE who gets higher throughout 
from a pico gets associated with a pico. In the next step 
called ABS rounding, the UE association decisions are used 
to obtain the ABS and non-ABS subframes. Indeed, this step 
produces a feasible A^'s and A/'^'s as we show later in our 
result. Finally, in the throughput computation step, each UE's 
available average airtime is scaled to fill-up the available sub- 
frames. The throughput of each UE can be computed from 
this. 

Performance guarantees: The worst case performance guar- 
antee of the output produced by Algorithm [2] depends on 
the number of iterations and step-size used for running Al- 
gorithm [T] prior to running Algorithm [2] This is shown by 
the following result. Let i?* be the throughput computed by 
Algorithm |2] and let be the optimal throughput. 

Proposition 3. Algorithm |2] produces feasible output to the 
problem OPT-ABS. Furthermore, for any given (5 > 0, there 
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Algorithm 2 Round RELAXED-ABS: Algorithm for Inte- 
ger Rounding of Output of Algorithm [T] 
1: User Association: For all u ^U, perform the following steps: 

1) Compute the throughput u gets from macro and pico in 
RELAXED-ABS solution as follows: 

rymacro macro ^ 

rypico pico,ABS I pico ^nA 

-^u ^li Vu ~r f'u Vu 5 

where, Xu^Vu i Vu^ are the out of Algorithm [l] 

2) If i?^^^^^ > i?^^^ UE-i^ associates with the macro, else 
it with pico. 

Define and compute U^, the set of UEs associated with macro- 
m after the UE association step. Similarly define and compute 
Up for every pico-p. 
2: ABS Rounding: Compute integral A^^'s and A*'s as follows: 

N:^ = Rndiv,^ (iVm) "im^M 
a; =Rndiv,^(ip) VpGP 

where Nm and Ap denote the output of Algorithm [ll 
3: Throughput computation: For each macro-m, for all G U^, 
the final value of x* , i?* are 



and Rl = 



(24) 



where Xm - J2ueu^ • 

For every pico, compute the ABS-utilization and non-ABS 



utilization Yp^ as 



.A ^nA .nA 



Next, for each pico-p, for all u G U^, the final values of 

Vu^^Vu^^^K are 

A* _ nA- _ K^jN^f-A;) 

P P 
7-)* pico.ABS A* I pico nA* /o/cn 

The system utility is computed as Util(R*) = Wu Ini?* . 



exists T large enough (but, polynomial in the problem param- 
eters and 1/5) and 7 satisfying ( [27] ) such that, if we apply 
Algorithm^to the output of Algorithm with this T, 7, then 



A. Cell Selection Bias for UE Association 

Our solution so far solves the coupled problem of optimiz- 
ing ABS sub-frames and UE association. However, UE asso- 
ciation in LTE HetNets need to be standard compliant. While 
standards on this are evolving, one proposed methodology by 
the LTE standard is a rule based on cell selection bias Q. 
Precisely, if be is the cell-selection bias of cell-c (which could 
be a macro or apico), then UE-^^ associates with cell such 
that, 

Ca = max[6c + RSRPu,c\ , (27) 

c 

where RSRPu,c is the received RSRP of cell-c at UE-^x. The 
choice of is a design issue. We wish to choose values of 
be so that UE association based on rule given by (27) leads 
to association decisions derived in the previous section by our 
optimization algorithm. 

Obtaining cell specific biases that precisely achieves a 
desired association may not always be feasible. Thus, we 
propose to compute biases so that the "association error" (as 
compared to optimal association) is minimized. We do this 
using the following steps. 

1) Since the relative bias between picos and macros matter, 
set all biases of macros to zero. 

2) Let Cp^rn be the set of UEs who have pico-p as the best 
candidate pico, and also, macro-m as the best candidate 
macro, i.e.. Cry rn = H Um- From the UEs in the set 



Cp^rn^ let ^p,m the total weight of UE's associated to 
pico-p under UE association produced by our algorithm 
in the previous section. Also from UEs in Cp^rn, as a 
function of bias 6, let Wp^rni^) be the total weight of 
UEs that would associate with pico-p if the bias of pico- 



p were set to b. In other words, Wp^rn = 
where 



Dp,m{b) — {it G Cp^m ' RSRPu^p -\- b > RSRPu,m} 9 

where RSRPu^m and RSRPu^p are received power (in 
dBm) of reference signal at UE-ii from best candidate 
macro and best candidate pico respectively. This step 
computes I^*^ and Wp,m(^) for every interfering pico- 
macro pair (j?, mO and every permissible bias value. 
3) For every pico-p, cell selection bias bp is set as 



Util{2{l + (5)R*) > Util{IC^^) 



Proof See (51. 



□ 



Remark 6. Proposition^shows that, for sufficiently large but 
polynomial number of iterations, the worst case approximation 
factor is close to 2. It is important to realize that, as with 
all NP-hard problems, this is simply a worst case result. 
Our evaluation with several real topologies suggest that the 
performance of our algorithm is typically within 90% of the 
optimal. Also, in practice, we recommend using the step-size 
and number of iterations given by ([22]). 



VIIL Computing Cell Selection Bias and ABS 
Patterns 

In this section, we describe two important computations that 
are relevant for realization of elCIC, namely. Cell Selection 
Bias (CSB) based UE association and converting ABS num- 
bers into ABS patterns. 



bp = arg mm 

h 



m^Xp 



(28) 



Thus, the bias values are chosen as the one that mini- 
mizes the mean square error of the association vector of 
number of UEs to different picos. 

Remark 7. ('Maximum and minimum bias constraint, j 
In many scenarios, operators that deploy picos may desire to 
have a maximum or minimum bias for a pico-p (say, bp^^nax 
and bp^rnin)' Por example, if a pico is deployed to fill a 
coverage hole or high demand area, then the bmin should be 
such that UEs around the coverage hole get associated with 
the pico. This can be handled using the following steps: 
1) Pirst run the joint UE-association and ABS- 
determination algorithm (Algorithm^and Algorithm^ 
by setting r^^^^^ = for all UEs that get associated 
with the pico even with minimum bias, and r^*^^ = 
for all UEs that do not get associated with the pico 
even with maximum bias. 
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2) Next execute the 3 steps for cell-bias determination 
described in this section but with the minor modification 
in ([28]) so that the argmin operation is restricted to 

B. Converting ABS numbers into ABS patterns 

In the previous sections, we have provided techniques to 
compute number of ABS subframes for every Nsf subframes 
(i.e., ABS subframes per ABS-period). In practice, we also 
need to specify the exact subframes in an ABS-period that are 
used as ABS subframes. The ABS number can be converted 
into a pattern as follows: 

1) Index the subframes in an ABS-period in a consistent 
manner across all macros and picos. 

2) Suppose a macro-m leaves out k out of Nsf subframes 
as ABS subframes. Then macro-m offers the first- k 
subframes as ABS subframes where the first- k relates 
to the indexing in the previous step. 

This simple scheme works provided all macros have the 
same set of permissible sub-frames if required (i.e., there is 
no restriction on certain macro that it cannot offer certain 
subframes as ABS). Notice that, since a pico can effectively 
use the least number of ABS offered by interfering macros, 
this scheme would naturally ensure a provably correct mapping 
between number of ABS and ABS-pattern. 

IX. Evaluation using RF Plan from a Real 
Network 

We evaluated our algorithms using RF plan from a real 
network deployment by a popular operator in Manhattan, 
New York City. The goal of our evaluation is four folds. 
First, to compare our proposed elCIC algorithm with other 
alternative schemes. Second, to understand the optimality gap 
of our algorithm because we have shown that optimizing 
elCIC parameters is NP-hard. Third, to understand the benefits 
offered by elCIC because operators are still debating whether 
the additional complexity of elCIC is worth the gains. Fourth, 
we show some preliminary results on how elCIC gains vary 
with pico transmit power and UE density. 

A. Evaluation Framework 

Topology: We used an operational LTE network deployment 
by a leading operator in New York City to generate signal 
propagation maps by plugging in the tower and terrain infor- 
mation along with drive-test data into a commercially available 
RF tool that is used by operators for cellular planing |1 1. In 
Figure [3] we show the propagation map of the part of the 
city that we used for evaluation, along with the macro-pico 
interference graph for nominal pico transmit power of 4W. 
The RF plan provides path loss estimates from actual macro 
location to different parts of the city. For the purpose of this 
study, we selected an area of around 8.9 km in the central 
business district of the city. This part of the city has a very 
high density of macro eNB's due to high volume of mobile 
data-traffic. The macro eNB's are shown in blue color with 
sectorized antennas and these eNB's are currently operational. 
While macro cells used in our evaluation are from the existing 
network, LTE pico cells are yet to be deployed in reality. Thus 
the pico locations were manually embedded into the network 
planning tool. We carefully chose 10 challenging locations for 
our picos: some are chosen with locations with poor macro 
signals, some pico locations are chosen with high density of 
interfering macros, some are chosen to coincide with traffic 
intensity hotspots, and one pico is also deeply embedded 




Pico^' ■ 




Fig. 3. Propagation map of the evaluated LTE network in New York City and 
the associated macro-pico interference graph. The blue eNB's are the macros 
and are currently operational. The grey eNB's are low power picos and we 
manually placed them in the tool. Pico- 10, Pico-3, Pico-5, Pico-9 have traffic 
hotspots around them. 

into a macro-cell. The picos are shown in red circle with 
omnidirectional antennas. 

All eNB's support 2x2 MIMO transmissions. UEs have 
2x2 MIMO MMSE-receivers that are also equipped with in- 
terference cancellation (IC) capability to cancel out broadcast 
signals from macro during pico downlink transmissions over 
ABS subframes. 

Important Cell Parameters: The macro eNB's have trans- 
mit power of 45 dBm (31 W). For the picos, we evaluated 
with 3 different settings of transmit powers: 36 dBm (4 W), 
30 dBm (1 W), and 27 dBm (500 mW). The bandwidth is 
10 MHz in the 700 MHz LTE band. The pico heights are 
chosen as typically 30 ft above the ground and macro heights 
are variable based on actual deployment and are typically 
much higher (more than 100 ft in many instances.). We choose 
Nsf = 40 so that we can obtain the number of ABS offerings 
in every 40 subframes. Also, we allow a maximum bias of 
15 dB for any pico because this is a typical restriction in 
current networks. 

Traffic: While the macros and the propagation map used in 
our evaluation is for a real network, we create synthetic UE 
locations for our evaluation because LTE pico deployments are 
still not very prevalent. This is done as follows. In the area 
under consideration, we chose a nominal UE density of around 
450 active UEs/sq-km (dense urban density). In addition, we 
created UE hotspots around Pico- 10, Pico-3, Pico-5, and Pico- 
9. The hotspots around Pico-3, Pico-5, Pico-9 have double the 
nominal UE density and the traffic hotspot around Pico- 10 
has 50% more UE density than nominal. We also performed 
evaluation by varying the UE density around the macro cells 
to 225 active UEs/sq-km (urban density) and 125 UEs/sq-km 
(sub-urban density) without altering the hotspot UE densities 
around the selected picos. As we discuss in Section |Xj in 
practice, network measurements would be available in terms 
of average traffic load and SINR distribution from which the 
UE locations can be sampled. 

Methodology: The radio network planning tool (RNP) 111 
and our elCIC implementation were used to generate the 
results as follows. 

1) RNP tool was used to generate signal propagation matrix 
in every pixel in the area of interest in New York City as shown 
in Figure [3] 

2) The RNP tool was then used to drop thousands of UEs 
in several locations based on the aforementioned UE density 
profiles. All UEs had unit weights (as in weighted proportional 
fair). 

3) Based on the signal propagation matrix and UE locations, 
we invoked the built-in simulation capability of RNP tool to 
generate the macro-pico interference graph and the following 
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donwlink SINR's for every UE: best-macro SINR, best-pico 
ABS SINR, best-pico non-ABS SINR. These SINR values 
were converted to physical layer rates using LTE look-up table. 
This step essentially produces the complete set of inputs for 
OPT-ABS problem as follows. 

4) Then this input-data was fed into our implementation of 
propose d elCI C and other comparative schemes described in 
Section [IX^ 

Thus we used RNP too to generate synthetic input that is 
representative of SINR and path-losses in a live network. 

B. Comparative elCIC Schemes 

The three schemes we compare are as follows. 

1) Proposed elCIC: This is the proposed algorithm de- 
veloped in this paper. Just to summarize, we first apply 
Algorithm [T] and the rounding scheme in Algori thm [2| a nd 
finally we use the technique described in Section VIII-A for 
obtaining CSB's. 

2) Fixed elCIC Pattern: Another option is to use a fixed or 
uniform elCIC pattern across the entire network. In iTO"], the 
authors have performed evaluation with fixed elCIC patterns. 
Also 1 12] considers fixed elCIC parameters. We also compare 
our proposed elCIC algorithm to the following four (ABS, 
CSB) combinations: (5/40, 5 dB), (10/40, 7.5 dB), (15/40, 
10 dB), (15/40, 15 dB). The fixed patterns represent the range 
of elCIC parameters considered in the literature. 

3) Local Optimal Heuristic: This is a local optimal based 
heuristic that is very easy to implement and is also amenable 
to distributed implementation. This scheme works as follows. 
First, each pico sets individual biases to maximize the total 
improvement (as compared to zero-bias) of physical layer rates 
(by considering the ABS rates) of all UEs within the coverage 
range of the pico. This step readily provides the set of UEs that 
associate with picos. In the next step, each macro m obtains 
the fraction (say, am) of UEs within its coverage range that 
associates with itself and then the macro offers \Nsf{l — am)] 
as ABS sub-frames. Each pico can only use minimum number 
of ABS sub-frames offered by its interfering macros. 

C. Results 

For our results, we consider all UEs in the coverage area 
of deployed picos and all macros that interfere with any of 
these picos. Clearly, these are the only UEs that are affected 
by elCIC or picos. 

Comparison with other schemes: In Figure |4] we compare 
our algorithm to different network wide fixed ABS settings. 
The interesting comparison is between our proposed elCIC 
algorithm and fixed (ABS, CSB) setting of (15/40, 15 dB) 
which corresponds to a CSB value of maximum possible 
15 dB; all other fixed elCIC schemes perform poorly. This 
fixed elCIC setting of (15/40, 15 dB) appears to perform better 
than our scheme for all UEs in the pico footprint area because 
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TABLE II 

Comparison of total log-utility (total of logarithm of UE 

THROUGHPUTS IN KBPS/MHZ) FOR DIFFERENT MACRO UE DENSITY AND 
PICO POWERS. DU, U, SU STAND FOR DENSE URBAN, URBAN, AND 
SUB-URBAN UE DENSITY. 
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2124.5 
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it associates all UEs in the pico footprint area to the pico; 
whereas, our scheme does not necessarily associate all UEs in 
the pico footprint area to the pico. However, the fixed elCIC 
scheme fails to account for the overall network performance 
as the macro UEs have to sacrifice far greater (compared to 
our scheme) throughput due to elCIC. Indeed, the throughput 
percentiles of all UEs in the system, for any fixed elCIC 
scheme, is reduced compared to our algorithm as we can see 
from the plot in the bottom panel of Fugure |4] For example, 
our proposed elCIC improves the 5^^, 10^^, 25^^ percentile 
throughput of Fixed-( 15/40, 15 dB) elCIC configuration by 
30 — 40%. For specific macros that do not interfere with 
any hotspot picos, this improvement is more than 50%. In 
typical deployments where many macros may have few or no 
pico neighbors in the macro-pico interference graph unlike our 
evaluation topology, the loss of the overall system performance 
could be more pronounced due to fixed elCIC configuration. 
Also, finding a good but fixed elCIC setting could also be 
challenging. Table |ll| also shows that the overall log-utilty of 
the system is better with our proposed scheme compared to 
the fixed elCIC schemes. 

In Figure [5j we compare the proposed the elC IC with the 
local optimal heuristic described in Section |IX-B[ Our scheme 



outperforms the local optimal heuristic by a margin of more 
than 80% for UEs in pico-footprint area; furthermore, the 
overall systems performance is better with our scheme as can 
be seen from the plot in the bottom panel of Figure [5] and 
Table |ll| However, the local optimal heuristic is very easy 
to implement and could be promising with additional minor 
changes. We leave this as a future topic of research. 

Optimality gap of our algorithm: Since the solution to 
RELAXED- ABS is an upper bound to the optimal solution 
of OPT-ABS, we obtain the optimality gap by comparing 
our final solution to that produced by RELAXED- ABS ( 
Algorithm [T]). We compute g, such that our algorithm is within 
100 X (1 — ^)% of the optimal, as follows. Suppose R^^^ 
and R^^^ be the UE-ii's throughputs produced by RELAXED- 
ABS and our complete algorithm respectively. Then, we say 
that the optimality-gap is a factor ^ < 1 if ln(i?^^^) > 
\n{Rlf^{l — g)). The smallest value of g that satisfies this 
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can easily be computed. In Table [Ill| we show for various 
settings of macro UE densities and pico transmission power 



Fig. 7. % -Improvement in 5*^,10*^,50*^ percentile of UE throughput 
with elCIC as opposed to no elCIC based pico deployment for different pico 
transmits powers (left-panel plot) and different macro cell densities (right- 
panel plot). 
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Benefits of elCIC: In the interference graph, we have 26 
macros and 10 picos. In typical deployments, there are going to 
be many more picos and macros. Thus, to understand the gains 
that even a few picos can offer, we show the following plots 
in Figure [6j CDF of throughputs of UEs in the pico coverage 
area, and CDF of throughput of UEs outside of pico coverage 
area. Thus, we wish to understand the gains of UE who 
could potentially associate with the picos, and the performance 
impact of UEs who do not have the option of associating 
with picos. The plot in the top panel of Figure [6] shows the 
throughput gains: (i) compared to no elCIC based scheme 
the gains are more than 200% for the far-edge UEs (say, 
2.5^^ percentile of the throughputs) and 40-55% for edge UEs 
(5^^ — 10^^ percentile of UE throughput), (ii) also, compared to 
no pico, the gains are even more dramatic and around 300% 
even for 5^^ percentile of the throughputs. The plot in the 
bottom panel of Figure [6] shows that the throughput gains 
(over no elCIC based pico deployment) of pico UEs do not 
come at an appreciable expense of macro UEs' throughput. In 
other words, though the macro eNB's have fewer subframes for 
transmissions (due to ABS offered to picos) using elCIC, this 
is compensated by the fact the macro UEs compete with fewer 
UEs (many UEs end up associating with picos under elCIC). 
Thus, there are great benefits of not only pico deployments, 
but also elCIC based pico deployments. 

elCIC gains witii power and load variation: To better 
understand the elCIC gains, in Figure [7] we compare the gains 
of elCIC using our algorithm with pico deployment without 
elCIC by varying the pico transmit powers and macro UE 
densities. In the top panel, we show the percentage throughput 
gain of elCIC scheme for different pico transmit powers for 
5^^, 10^^, 50^^ percentile of UE throughputs. It can be seen 
that, it is the edge UEs who really gain with elCIC; indeed, 
this gain could even come at the expense of UEs close to the 
pico (as can be seen with 1 W pico power scenario) who do 
not gain much due to elCIC. The edge gain is also a direct 
consequence of our choice of log-utility function as system 
utility. In the bottom panel of Figure [7] we show the gains for 
different macro UE density. It can be see that, higher macro 
UE density results in higher gain due to elCIC. Intuitively 
speaking, more the UEs that have the choice of associating 



with picos, larger are the elCIC gains from our algorithm. This 
suggests the usefulness of our scheme for practical scenarios 
with large number of picos and very high density in the traffic 
hot-spot areas. 

Optimal Parameters: In Table IV we show the optimal 
ABS received by each pico and the associated bias obtained 
using our algorithm. There are a couple of interesting obser- 
vations. First, the ABS offered to picos not only depends on 
traffic load but also depends on the number of interferers. For 
example, Pico-5 received 4 out 40 subframes for ABS though 
it has a hotspot around it, however, Pico-1 received 6 out of 
40 subframes for ABS without any hotspot around it. This is 
because, Pico-5 has more neighbors in the macro-pico interfer- 
ence graph. Second, Pico-2 which is embedded into a macro, 
also receives 4 ABS subframes and serves as enhancing in-cell 
throughput. Thus picos can go beyond improving throughput 
in edges if elCIC parameters are configured in a suitable 
manner. This also shows that there could be considerable 
variation in optimal ABS and CSB settings. This explains the 
poor performance of network wide fixed elCIC schemes. 

X. SON AND ElCIC: Challenges and Discussion 

A key aspect of LTE networks is its Self Optimized Net- 
working (SON) capability. Thus, it is imperative to establish 
a SON based approach to elCIC parameter configuration of 
an LTE network. The main algorithmic computations of SON 
may be implemented in an centralized or a distributed manner. 
In the centralized computation, the intelligence is concentrated 
at the Operations Support System (OSS) layer of the network, 
while in the distributed computation the computation happens 
in the RAN or eNB. The main benefit of a centralized approach 
over a distributed approach is twofold: a centralized solution 
in OSS is capable of working across base stations from 
different vendors as is typically the case, and well-engineered 
centralized solutions do not suffer from convergence issues of 
distributed schemes (due to asynchrony and message latency). 
Indeed, realizing these benefits, some operators have already 
started deploying centralized SON for their cellular networks. 
Nevertheless, both centralized and distributed approaches have 
their merits and demerits depending on the use-case. Also, it is 
widely accepted that, even if the key algorithmic computations 
happen centrally in OSS, an overall hybrid architecture (where 
most heavy-duty computations happen centrally in OSS with 
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distributed monitoring assistance from RAN) is best suited for 
complicated SON use-cases such as elCIC whereas complete 
distributed approach is suited for simple use-cases like cell- 
neighbor detection. We next discuss our prototype hybrid SON 
for elCIC. 

A. Hybrid SON Architecture 

The architecture of our prototype is shown in Figure [8] 
Apart from an operational wireless network, the architecture 
has two major component-blocks: a network planning tool and 
an engine for computing optimal elCIC configuration. 

In our prototype the main optimization task is executed cen- 
trally at the OSS level, but input distributions are provided by 
the RAN and further estimated/processed in the OSS. From the 
perspective of the optimization algorithm, centralization offers 
the best possible globally optimal solution given accurate and 
timely data inputs. We also assume that other affiliated SON 
procedures, such as Automatic Neighbor Relations (ANR) are 
executed in the RAN and their results reported to the OSS. 
After the execution of the OSS optimization algorithms, the 
optimal parameters propagate southbound towards the RAN 
elements. 

Monitoring in the RAN: As part of the Operations Admin- 
istration and Maintenance (OA&M) interfaces. Performance 
Management (PM) data are reported to all OSS applications, 
SON-applications including. In general, we can have periodic 
reporting or event-based reporting, an implementation choice 
that trades latency and accuracy. Irrespective of the imple- 
mentation, elCIC requires from the RAN, path loss statistics, 
traffic load statistics and SINR statistics. Our prototype im- 
plementation can flexibly accommodate the case of missing 
data i.e. incomplete statistics, that is common during the 
planning phase of the network or in the case where the required 
inputs are not readily available by another vendor's RAN 
implementation. In these cases we can easily replace actual 
network data with synthetic data generated by a radio network 
planning tool. 

The role of radio network planning (RNP) tools: Var- 
ious Databases are used to import information necessary for 
performing the radio network planning. Inventory information 
that provides network topology, drive-tests that calibrate path 
loss models as well as performance measurement data that 
determine the shape and value of traffic intensity polygons in 
the area of interest, are the most important information sources 
aggregated in the RNP tool. RNP can use this information to 
generate synthetic input data (using built-in simulation capabil- 
ities of the tool) for our elCIC algorithm. For our purpose, we 
used a planning tool 1 1 1 that uses the traffic map, propagation 
map, and eNodeB locations to generate multiple snapshots of 
UE locations. This input data is saved into a database. The 
elCIC computation engine implements our proposed algorithm 
to compute optimal elCIC parameters. The role of the radio 
network planning tool as a prior distribution generator is 
quite important in the planning as well as the initialization 
procedures of a network element when measurements are 
unavailable or too sparse. 

B. Computational Flow 

Average Load Based Input: The traffic data that is avail- 
able from the real network comprises of average traffic load 
information over a period of interest. Our method can be 
used for the purpose of obtaining the correct ABS and CSB 
configuration can be used as follows. 
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Fig. 8. Prototype Architecture. 

We consider a scenario where, for each cell, average traffic 
load and SINR distribution (in different subframes including 
ABS subframes) are reported periodically. Indeed, such report- 
ing is common in many real network deployments 1 151. We 
adapt our solution as follows: 

1) Generation of sample UE location: Based on the traffic 
map and SINR distribution, multiple system wide sam- 
ple UE-location snapshots are generated (commercial 
network planning tools usually have this capability). 
Each sample UE-locations are translated into downlink 
PHY-layer rate between UE and Macro, UE and Pico 
with and without ABS. Note that, the exact location of 
a sample UE is of little relevance, rather, the SINR to 
the nearest macro and pico along with the RSRP values 
are of relevance here. 

2) Optimization for each snapshot: Based on the UE- 
location snapshots, for each such sample snapshot, our 
algorithm is run for solving OPT-ABS. This generates 
ABS and CSB configuration for each sample. 

3) Monte-Carlo averaging: Once ABS and CSB compu- 
tation is performed for sufficient sample UE location 
snapshots, the results are averaged over all samples. 

Time-scale and dynamic elCIC: There are two important 
considerations for the reconfiguration frequency of elCIC 
parameters. Firstly, since only the traffic distribution statistics 
can be obtained from the network, it is imperative that elCIC 
computations happen at the same time-scale at which the traf- 
fic distribution can be estimated accurately; otherwise, elCIC 
changes at a faster time- scale may not provide appreciable 
gains while causing unnecessary reconfiguration overhead. 
Secondly, since elCIC reconfiguration involves a cluster of 
macros and picos, it takes a few minutes to have new set of 
traffic information from all the cells ifTSl . Therefore, the time 
scales of changing elCIC parameters for all practical purposes 
are in the order of few minutes (typically, 5 — 15 minutes). 
Thus, elCIC reconfiguration ought to happen whenever (i) 
traffic load changes significantly in some cells, or ( ii) or when 
a maximum duration elapses since the last reconfiguration. In 
addition, if the estimated improvement upon new elCIC re- 
computation is small, then the previous elCIC configuration 
can be maintained. 

C. Fully Distributed SON Architecture 

Another implementation option is to distribute elCC compu- 
tations at the Network Element (eNB) level and evidently what 
is traded off with this approach, as compared to centralized, is 
signaling and communication latency. In LTE, the X2 interface 
can be used to interconnect eNBs and this interface has been 
the subject of extensive standardization when it comes to 
interference management. The main challenge in distributed 
approach comes from X2 latency and asynchrony. 
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Using proprietary messages over X2, our solution can be 
adapted for in-network computation where the macro and 
picos use their local computational resources to obtain a 
desired solution. We provide a broad outline in the following. 
The key to achieving this is our dual based implementation 
of ABS-RELAX. Suppose the optimization task has to be 
performed periodically or over a time-window of interest. 
Then, a distributed implementation involves the following 
high-level steps, (i) Macros and picos generate UE samples 
in the cells based on traffic/SINR distribution of UEs in the 
time- window of interest; the sample data is exchanged between 
neighboring macro-pico pairs. Algorithm [T] is run in a 
distributed manner, where, in each iteration, macros and picos 
exchange dual variables jj^p^rn and A^'s for relevant UEs. This 
exchange allows macros and picos to update primal variables 
(for themselves and also candidate UEs) locally, (in) Finally, 
note that the rounding step and UE association in Algorithm |2] 
can be performed locally. If UE association is implemented via 
cell- selection bias, the scheme described in Section [VIII-A| can 
be easily carried out locally at each pico (this may require 
collecting UE association vector by message exchange with 
neighboring macros). 

XI. Concluding Remarks 

In this work, we have developed algorithms for optimal 
configuration of elCIC parameters based on actual network 
topology, propagation data, traffic load etc. Our results using 
a real RF plan demonstrates the huge gains that can be had 
using such a joint optimization of ABS and UE-association 
based on real network data. The broader implication of our 
work is that, to get the best out of wireless networks, networks 
must be optimized based on real network data. 
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